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Abstract. The Universe on large scales is well described by the ACDM cosmological model. 
There however remain some heavy clouds on our global understanding, especially on galaxy 
scales, which we review here. While some of these clouds might perhaps disappear through 
small compensatory adjustments of the model, such as changing the mass of the dark matter 
particles or accounting better for baryonic physics, others should rather be taken as strong 
indications that the physics of the dark sector is, at the very least, much richer and complex 
than currently assumed, and that our understanding of gravity and dynamics might also be at 
play. For instance, the empirically well-tested MOND phenomenology in galaxies, whatever its 
final explanation, should be understood in any model of galaxy formation and dynamics. Current 
alternatives to ACDM however bring with them many unsolved questions and challenges. 



1. Introduction 

Assuming General Relativity to be valid on all scales, data ranging from the Cosmic Microwave 
Background to individual galaxies point towards a Universe dominated by dark energy and dark 
matter, the nature of these being most certainly amongst the deepest problems of modern 
physics. While dark energy is well-represented by a cosmological constant A in Einstein's 
field equations, the currently preferred dark matter candidate is a collection of stable, neutral, 
elementary particles that condensed from the thermal bath of the early Universe, and which are 
known as 'cold dark matter' (CDM) particles (see, e.g., Frenk & White 2012, Strigari 2012 for 
recent reviews). On galaxy scales, however, predictions of this standard ACDM cosmological 
model, although plagued by the enormous complications of baryonic astrophysics, are difficult 
to reconcile with observations. We hereafter review these challenges for the ACDM model, and 
point out that some of them hint at a richer and more complex physics of the dark sector 
than currently assumed. In particular, many observed scaling relations involve the ubiquitous 
appearance of an acceleration constant ao ~ A 1 / 2 ~ 10~ 10 m/s 2 , whose origin is a deep mystery 
in the standard context. Surprisingly, most of these scaling relations can be summarized by the 
empirical formula of Milgrom (1983). The success of this formula means that the gravitational 
field in galaxies mimicks, for whatever reason, an effectively modified force law on galaxy scales, 
known as Modified Newtonian Dynamics (MOND). This formula however fails to account for 
dynamics and lensing of galaxy clusters (e.g., Clowe et al. 2006, Angus et al. 2007), meaning 
that if the formula would be reflecting a true modification of gravity as an alternative to galactic 
dark matter, it should rely at cluster scales on residual missing mass, which could be in baryonic 



or non-baryonic form (Milgrom 2008, Angus et al. 2007), or the formula should be extended 
(Zhao & Famaey 2012). It is also unclear how the angular power spectrum of the Cosmic 
Microwave Background (CMB) could be accounted for without resorting to some form of non- 
baryonic dark matter (e.g., Slosar et al. 2005, Angus 2009). Nevertheless, the main motivation 
for studying alternatives to ACDM based on the MOND phenomenology is not necessarily to 
get rid of any form of 'dark matter', but rather to explain why the observed gravitational field 
in galaxies is apparently mimicking a universal force law generated by the baryons alone. The 
simplest explanation is of course a priori not that dark matter particles arrange themselves 
(by a hitherto unknown physical mechanism) in order to mimick a fake force law, but rather 
that the force law itself is modified. However, at a fundamental level, covariant theories of 
modified gravity often have to include new fields in the dark sector to reproduce this effective 
force law (fields with an energy density nevertheless subdominant to the baryonic one, and a role 
completely different from that of CDM: they would mediate the force rather than sourcing it), or 
even introduce what could be called a 'dark matter medium' (with an energy density outweighing 
the baryonic one) exhibiting a kind of new fundamental interaction with baryons: this makes 
the confrontation between MOND and dark matter much less clear than often believed, since 
the former implies a more complex structure of the dark sector than the currently assumed 
CDM, but does not necessarily imply the absence of a dark sector. In MOND, the new fields 
responsible for the effects usually attributed to CDM would in fact be somewhat closer to dark 
energy fields than to CDM. It would of course be even more exciting if one would manage to 
find a physical connection between these putative new fields and the actual dark energy sector. 

We herefater list a (non-exhaustive) series of problems for ACDM (Sect. 2), then show how 
theories based on the MOND phenomenology might address a lot of them (Sect. 3), before listing 
a (non-exhaustive) series of covariant theories currently proposed in this vein (Sect. 4), together 
with their own questionmarks and internal problems. The reading of this short proceeding can 
be completed by the reading of the exhaustive review paper recently published in Living Reviews 
in Relativity (Famaey k McGaugh 2012). 

2. Current problems for our understanding of galaxy properties in ACDM 

2.1. The nature and distribution of dwarf galaxies 
2.1. a. The missing satellites challenge. ACDM simulations predict vast numbers of subhaloes 
that are satellites to the main halo hosting a galaxy like the Milky Way. One would naively 
expect each of these subhaloes to host their own miniature satellite galaxy. The Local 
Group looks nothing like this prediction, having only a small handful of dwarfs around each 
of the giant galaxies. Taking into account stellar feedback and heating processes (mainly 
at re-ionisation) in the galaxy formation process, the predicted number of faint satellites 
around a Milky Way-like galaxy is ~ 100 to 600. Since the majority of the 24 known 
satellites of the Milky Way have been largely discovered with the Sloan Digital Sky Survey 
(SDSS), and since this survey covered only one fifth of the sky, it has been argued that 
the problem was solved. However, models that successfully explain the properties of the 
Milky Way dwarf satellites predict unobserved dwarfs beyond the virial radius (Bovill &; 
Ricotti 2011a). Moreover, the models tend to produce an overabundance of bright dwarf 
satellites (Ly > 1O 4 L0) with respect to observations (Bovill & Ricotti 2011b). A rather 
discomforting way-out is to simply state that the Milky Way must be a statistical outlier, but 
this is actually contradicted by the study of Strigari & Wechsler (2012) on the abundance of 
bright satellites around Milky Way-like galaxies in the SDSS survey. A somewhat related but 
nevertheless distinct problem is that simulations predict that the most massive subhaloes of 
the Milky Way are too dense to host any of the observed bright satellite galaxies (Boylan- 
Kolchin et al. 2011, 2012). This is the 'too big to fail' aspect of the missing satellites 
challenge. A potential solution might be warm rather than cold dark matter. 



2.1.b. The satellites phase-space correlation challenge (disks of satellites). The 

distribution of dark subhalos around Milky Way-sized halos is also predicted by ACDM to 
be roughly isotropic. However, the Milky Way satellites are currently observed to be highly 
correlated in phase-space: they lie within a seemingly rotationally supported, relatively thin 
disk (Fig. 1, see Kroupa et al. 2010). This differs from the prediction for subhalos in both 
coordinate and momentum space. Since the SDSS survey covered only one fifth of the sky, 
it will of course be most interesting to see whether ongoing surveys such as Pan-STARRS 
will confirm this state of affairs. Whether or not such a satellite phase-space correlation 
would be unique to the Milky Way should also be carefully checked: the Milky Way could 
be a statistical outlier, but if the ACDM model is a realistic description of nature, then the 
average satellite configurations in external galaxies should be only moderately flattened, 
and generally non-rotating. In this respect, the recent findings of Ibata et al. (2013) that 
half of the satellites of M31 define an extremely thin and extended rotating structure makes 
this phase-space correlation challenge even more severe (see Sect. 3 for a possible solution). 



Figure 1. Spatial distribution of 
the Milky Way (MW) satellites (from 
Kroupa et al. 2010). Both axes in 
kpc. The 11 classical satellites are shown 
as large (yellow) circles, the 13 recently 
discovered satellites are represented by 
the smaller (green) dots, Pisces I and II 
being the two southern dots. The two 
open squares near the MW are Seg 1 and 
2. The obscuration-region from the MW 
disk is given by the horizontal gray area. 

2. 1. e. The density-morphology relation of dwarf ellipticals. More dwarf elliptical galaxies 

are observed in denser environments (e.g. Kroupa et al. 2010). This relation, observed in 
the field, in galaxy groups and in galaxy clusters, is not yet understood. 

2.2. The dynamical friction and stabilization from dark matter 

2. 2. a. The angular momentum challenge. Both the merger history of galaxy disks in a 

hierarchical formation scenario and the associated transfer of angular momentum from the 
baryonic disk to the dark halo cause the specific angular momentum of the baryons to end 
up being much too small in simulated disks. These simulated disks in turn end up much 
smaller than the observed ones, and elliptical systems end up too concentrated too. To 
address this challenge within the standard paradigm, one needs to form disks through late- 
time quiescent gas accretion from large-scale filaments, with less late-time mergers than 
currently predicted in ACDM. 

2.2.b. The pure disk challenge. Related to the previous problem, large bulgeless thin disk 
galaxies are very difficult to produce in simulations. This is because major mergers typically 
create bulges. This is fundamental to the hierarchical nature of galaxy formation in ACDM, 
wherein galaxies are built by the merger of smaller galaxies. Bulgeless galaxies should 
therefore represent the rare, quiescent tail of a distribution of merger histories for galaxies 
of the Local Volume. However, these bulgeless disk galaxies are quite common, representing 
more than half of large galaxies (with V c > 150 km/s) in the Local Volume. 
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2.2. c. The stability challenge. Quasi-spherical CDM halos stabilize very low surface density 

disks against the formation of bars and spirals, due to a lack of disk self-gravity. The 
observation (McGaugh et al. 1995) of Low Surface Brightness (LSB) disk galaxies with 
prominent bars and spirals is thus challenging in the absence of a significant disk component 
of dark matter. Moreover, in the absence of such a disky dark matter component, the lack of 
disk self-gravity prevents the creation of observed large razor thin LSB disks. This suggests 
more disk self-gravity than meets the eye (McGaugh & de Blok 1998). In the standard 
context, this would tend to point towards an additional disk-specific DM component, either 
a CDM-one linked to in-plane accretion of satellites, or a baryonic one, for instance in the 
form of molecular gas or due to a bottom-heavy IMF. 

2.3. Structure formation not fast enough in ACDM? 

2.3. a. The bulk flow challenge. In the ACDM model, peculiar velocities of galaxy clusters are 

predicted to be of order ~ 200 km/s. These velocities can actually be measured through 
the fluctuations in the CMB generated by the scattering of the CMB photons by the hot 
X-ray-emitting gas inside galaxy clusters. These observations yield an observed coherent 
bulk flow of order ~ 1000 km/s on scales out to at least 400 Mpc (Kashlinsky et al. 2012). 
This bulk flow challenge also appears in galaxy studies (Watkins et al. 2009). A related 
problem is the collision velocity of the merging bullet cluster 1E0657-56 at z = 0.3, which is 
larger than 3100 km/s, too high a value to be accounted for by ACDM. These observations 
seem to indicate that the attractive force is enhanced compared to what ACDM predicts. 
2.3. b. The high-z clusters challenge. The observation of a single very massive cluster at high 
redshift can falsify ACDM. In this respect the existence of galaxy clusters like El Gordo with 
a mass of ~ 2 x 10 15 M & at z = 0.87 (Menanteau et al. 2012) and XMMU J2235. 3-2557 with 
a mass of of ~ 4 x 10 14 M Q at z = 1.4 is surprising. Though not sufficient to rule out the 
model, these clusters certainly push the envelope: there is only about a 50:50 chance that 
a cluster as massive as El Gordo exists at such high redshift. This has been the empirical 
experience: for every generation of surveys, structures are found further away, and appear 
sooner than expected. 

2.3.c. The Local Void challenge. The 562 known galaxies at distances smaller than 8 Mpc from 
the center of the Local Group define the 'Local Volume'. Within this volume, the region 
known as the 'Local Void' hosts only 3 galaxies. This is much less than the expected ~ 20 
galaxies for a typical similar void in the ACDM model (Peebles & Nusser 2010). Moreover, 
in the Local Volume, large luminous galaxies are 6 times more abundant than predicted in 
the underdense regions. This could just mean that the Local Volume is a statistical outlier, 
but it could also point, in line with the two previous problems, towards more rapid structure 
formation in the early Universe, allowing sparse regions to form large galaxies cleaning their 
environment in a shorter time, thus making the voids emptier and the galaxies larger at 
early times. 

2.4- Systematics in the dark-to-baryonic matter ratios and the ubiquitous ao-scale 
2.4-a. The cusp-core challenge. A long-standing problem of ACDM is the fact that the 
numerical simulations of the collapse of dark matter halos lead to a 'cuspy' density 
distribution that rises steeply as a decreasing function of radius, while rotation curves 
of external galaxies imply nearly constant density cores in the central parts. The state-of- 
the-art solution to this problem is to enforce strong supernovae outflows that move large 
amounts of low-angular-momentum gas from the central parts and that pull on the central 
dark matter concentration to create a core, but this is still a relatively fine-tuned process, 
which fails to account for cored profiles in the smallest galaxies, and fails to produce their 
observed baryon fractions 



2.4-b. The missing baryons challenge(s): baryon fraction defined by the acceleration 
in units of do- Constraints from the CMB imply f2 m = 0.27 and f2& = 0.046. However, our 
inventory of known baryons in the local Universe comes up short of the total. For example, 
Bell et al. (2003) estimate that the sum of cold gas and stars is only ~ 5% of fV While 
it now seems that many of the missing baryons are in the form of highly ionized gas in the 
warm-hot intergalactic medium (WHIM), we are still unable to give a confident account of 
where all the baryons reside. But there is another missing baryons challenge within each 
CDM halo: one would indeed naively expect each halo to have the same baryon fraction as 
the Universe, /& = Qb/Q m = 0.17. On the scale of clusters of galaxies, this is approximately 
true (but still systematically low), but for galaxies, observations depart from this in a 
systematic way which remains unexplained in the standard context. In Fig. 2 (McGaugh et 
al. 2010), the ratio of the detected baryon fraction over the cosmological one, fy, is plotted as 
a function of the potential well of the systems. There is a clear correlation: at large radii, the 
baryon fraction is actually directly equivalent to the acceleration in units of ao = 10~ 10 m/s 2 . 
If we adopt a rough relation M500 — 1-5 x 10 5 M & x VJ?(km/s) -3 , we get that the acceleration 
at i?500i and thus the baryon fraction, is M^/M^qq = asoo/«o — 4 x 10 -4 x V^.(km/s) . 
Divided by the cosmological baryon fraction, this explains the trend for fd = M;,/(0.17M5oo) 
with potential (<3? = V^) in Fig. 2. This missing baryons challenge is actually closely related 



Figure 2. The fraction of de- 
tected over expected baryons fd = 
(Mb/Msoo)/ ft, as a function of po- 
tential well and mass. Each point 
represents many objects. Orange 
squares are Local Group dwarf 
satellites, light blue circles gas- 
dominated disk galaxies, dark blue 
circles star-dominated spiral galax- 
ies, grey triangles galaxy clusters. 

2.4-c The baryonic Tully Fisher relation (with Gao defining zero point). The baryonic 
Tully-Fisher relation is an observed relation between the asymptotic circular velocity and 
baryonic mass in galaxies (Fig. 3). Recent work extending the relation to low mass, typically 
low surface brightness and gas rich galaxies, extends the dynamic range of the relation to 
five decades in baryonic mass. Over this range, the relation has remarkably little intrinsic 
scatter (consistent with zero given the observational errors) and is well described as a power 
law, or equivalently, as a straight line in log-log space logM& = alogVy — log/3, with slope 
a = 4. This slope is consistent with a constant acceleration scale ao ~ Vjr / (GMf,) such that 
the normalization constant j3 = Ga^. All rotationally-supported galaxies are observed to 
follow the relation, irrespective of their size, surface brightness, or gas fraction. 

2.4-d Tidal dwarf galaxies. Even young tidal dwarf galaxies (TDGs) formed in the collision of 
larger galaxies appear to obey the baryonic Tully-Fisher relation (Gentile et al. 2007). This 
is surprising because these galactic collisions should be very effective at segregating dark 
and baryonic matter: the rotating gas disks of galaxies that provide the fodder for tidal 
tails and the TDGs that form within them initially have nearly circular, coplanar orbits. In 
contrast, the dark matter particles are on predominantly radial orbits in a quasi-spherical 
distribution. This difference in phase space leads to tidal tails that themselves contain 
very little dark matter. When TDGs form from tidal debris, they should be largely devoid 



to the baryonic Tully-Fisher relation. 
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of dark matter. Nevertheless, TDGs do appear to contain enough dark matter to obey 
the same baryonic Tully-Fisher relation as galaxies that have not suffered this baryon-dark 
matter segregation process. 




Figure 3. The Baryonic Tully- 
Fisher relation for galaxies with 
well-measured asymptotic veloci- 
ties Vf. The dark blue points have 
M* > M g , while the light blue 
points have M* < M g and are gen- 
erally less precise in velocity, but 
more accurate in terms of mass. 
The dotted line has the observed 
slope 4, the dashed line has slope 
3, the difference between these two 
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2.4- e. oo as a transition acceleration. The mass discrepancy in galaxies always appears 
(transition from baryon dominance to dark matter dominance) when V^/R ~ oq, yielding 
a clear mass-discrepancy acceleration relation (Fig. 4). This, again, is the case for every 
single rotationally supported system irrespective of its formation mechanism and history. 
For High Surface Brightness (HSB) galaxies, where there exists two distinct regions where 
V^/R > clq in the inner parts and V^/R < clq in the outer parts, locally measured mass-to- 
light ratios show no indication of hidden mass in the inner parts, but rise beyond the radius 
where V^/R « ao- It can never be emphasized enough that the role played by ao in the 
zero-point of the baryonic Tully-Fisher relation and this role as a transition acceleration 
have strictly no intrinsic link with each other, they are fully independent of each other 
within the CDM paradigm. There is nothing built in ACDM that stipulates that these two 
relations (the existence of a transition acceleration and the baryonic Tully-Fisher relation) 
should exist at all, and of course nothing that implies that these should exhibit an identical 
acceleration scale. 

2.4-f- clq/G as a transition central surface density. The acceleration ao defines the transition 
from LSB galaxies to HSB galaxies: in LSB galaxies whose central surface density is much 
smaller than some critical value of order S+ = clq/G, DM dominates everywhere, and the 
magnitude of the mass discrepancy is given by the inverse of the acceleration in units of 
ao- On the other hand, baryons dominate in the inner parts of HSB galaxies whose central 
surface density is higher than £+. There is an anti-correlation between the baryonic surface 
density and the fractional contribution of DM to the rotation curve (Fig. 5). To explain 
this, there must be a strong fine-tuning between dark and baryonic surface densities, a sort 
of repulsion between them, a repulsion which is however contradicted by the correlations 
between baryonic and dark matter bumps and wiggles in rotation curves (known as 'Renzo's 
rule'). The shapes of rotation curves also depend on surface density: HSB galaxies have 
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Figure 4. The mass discrepancy is 
defined as (V/Vb) 2 where V is the 
observed velocity and Vj, is the ve- 
locity attributable to baryonic mat- 
ter. Many hundreds of individ- 
ual resolved measurements along 
the rotation curves of ~ 100 spiral 
galaxies are plotted. The mass dis- 
crepancy is plotted as a function of 
radius (top panel), centripetal ac- 
celeration a = V 2 /r (middle panel), 
and Newtonian acceleration = 
V b 2 /r (bottom panel). 



Figure 5. The fractional contri- 
bution of baryons (Vb/V p , top) and 
dark matter (Vdm/Vp, bottom) to 
the rotation velocity V p at the ra- 
dius Rp where the contribution of 
the baryons peaks, plotted as a 
function of the characteristic sur- 
face density (E 6 = 0.75M fe /i?p. 
As the baryonic surface density in- 
creases, the fractional contribution 
of dark matter to the total gravi- 
tating mass decreases. This fine- 
tuning persists for any choice of 
stellar mass-to-light ratio. 



rotation curves that rise steeply then become flat, or even fall somewhat to the not-yet- 
reached asymptotic flat velocity, while LSB galaxies have rotation curves that rise slowly 
to the asymptotic flat velocity. Finally the total (baryons+DM) acceleration declines with 

1 /2 

the mean baryonic surface density of galaxies, in the form a oc S fe (Fig. 6). 
2.4-g- Features in the baryonic distribution imply features in the rotation curve. This 
is known as 'Renzo's rule' (Sancisi 2004, see also Swaters et al. 2012). While the effect 
of non-axisymmetric motions should be investigated in more detail, this is currently most 
easily interpreted by either a heavy 'dark baryonic' component (e.g., in the form of molecular 
H2) scaling with HI in galaxy disks, or by a theory where the baryons effectively act as the 
main source of gravity. 

2.4-h. ao/G as a critical mean surface density for stability. Disks with mean baryonic 
surface density (E) > £j = ao/G are extremely rare (Fig. 7) and unstable. If stability is 
provided by a DM halo, there is no reason for this stability threshold to be related to the 



same acceleration scale that appears elsewhere (so this is again an independent occurence of 
the ao-scale). Dense halos predicted by ACDM could (and should) in principle host higher 
surface brightness disks than are observed. 




Figure 6. The characteristic cen- 
tripetal acceleration a p = V p 2 / R p 
: in units of a® plotted against the 
: characteristic baryonic surface den- 
sity. The data do show a correla- 

1/2 

tion {a p oc Ti b ), clearly indicating 

ti 1 a dynamical role for the baryons, 

io~ 3 ict 1 io° io 1 whilst CDM dominance should nor- 
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Figure 7. The characteristic sur- 
face density of baryons vs. their 
dynamical scale length R p (Left), 
and disk-only central surface den- 
sity vs. disk scale-length (Right). 
High characteristic surface densities 
at low R p in the left panel are typ- 
ical of bulge-dominated galaxies, 
and disks with characteristic sur- 
face density > ao/G are extremely 
rare. 



3. Solving problems: Milgrom's law and MOND 

As we have seen, the ao acceleration-scale is ubiquitous in galaxies. It defines the zero point of 
the Tully-Fisher relation and the amplitude of the mass-discrepancy in the weak-gravity limit, it 
appears as a transition-acceleration above which no DM is needed, and below which DM appears, 
and it defines the maximum surface density of pure disks. It would obviously be interesting if 
these observational occurences could be summarized and empirically unified in some way: such 
a unification is actually feasible through the rather old idea of Milgrom (1983) that the total 
gravity g is related to the baryon-generated one gN in a way analogous to Coulomb's law in a 
dieletric medium: 

AM — Jg = gN, (1) 

where the interpolating function 

/j,(x) — > 1 for x 1 and fi(x) — >• x for x <C 1. (2) 

This law is at the basis of the Modified Newtonian Dynamics (MOND) paradigm. It does not 
simply make rotation curves flat, but predicted all the challenging observations of Sect. 2.4 
in a quick and efficient way (see Sect. 5.2 of Famaey &i McGaugh 2012). This formula was 
not constructed to explain already known observational facts, but predicted almost all of 



them a priori. Its salient properties can be elegantly summarized by simply stating that the 
underlying theory, whatever its deep nature, has to become scale-invariant under transformations 
(t,r) — > X(t,r) in the regime of small acceleration a <C oo (see, e.g., Milgrom 2012). 

Addressing all the other challenges, from Sect. 2.1 to 2.3, would need a full-fledged theory 
predicting Milgrom's law and the associated MOND phenomenology. As we shall see in the next 
section, such a full-fledged theory does not exist at the moment. One can however speculate on 
how these problems could be addressed by such a theory. For instance, if the fields generating 
the MOND behavior on galaxy scales have an energy density subdominant to the baryonic one, 
effects such as dynamical friction with DM would disappear. Galaxy disks would then easily 
get more extended thanks to the absence of angular momentum transfer to the dark halo, but 
mergers would also take longer, and be subject to multiple passages (Tiret & Combes 2008). 
While MOND would naturally provide added stability compared to Newtonian dynamics without 
DM, it does not over-stabilize disks: features like bars and spiral arms are a natural result of disk 
self-gravity, and numerical simulations based on MOND appear to do a good job at reproducing 
the range of observed morphologies of spiral galaxies (Tiret & Combes 2007, 2008). Concerning 
structure formation, the MOND paradigm has a priori no cosmology, providing analogs for 
neither the Friedmann equation nor the Robertson- Walker metric. For these, one must appeal 
to specific hypotheses for the covariant parent theory of MOND (Sect. 4), which is far from 
unique. What is more, none of the candidates is theoretically satisfactory, as none emerges 
from first principles. It is thus not clear whether a compelling candidate MOND cosmology 
will ever emerge. One can however speculate on what would be a reasonable scenario, based on 
simulations using MOND in the weak-field regime at their basis (e.g., Sanders 1998, Angus & 
Diaferio 2011). In the absence of CDM, and with a MOND-like force driving structure formation, 
after initially lagging behind, structure formation would rapidly speed up compared to ACDM 
once the influence of radiation declines and perturbations begin to enter the MOND regime. 
Large galaxies could form by z « 10 (Sanders 2008) and clusters by z 2, considerably earlier 
than in ACDM. By z = 0, the voids would become more empty than in ACDM, but otherwise 
MOND simulations (of collisionless particles, which is of course not the best representation of 
the baryon fluid) show the same qualitative features of the cosmic web (Llinares et al. 2008). 
The main difference is in the timing of when structures of a given mass would appear (assembling 
a large mass early being easier in MOND). This means that MOND is promising in addressing 
the high-z clusters challenge and Local Void challenge, as well as the bulk flow challenge and 
high collisional velocity of the bullet cluster (Llinares et al. 2009), due to the much larger than 
Newtonian MOND force in the structure formation context. What is more, it could allow large 
massive galaxies to form early (z ~ 10) from monolithic dissipationless collapse (Sanders 2008), 
with well-defined relationships between the mass, radius and velocity dispersion. Consequently, 
there would be less mergers than in ACDM at intermediate redshifts, but since these mergers 
would last longer, the number of interacting galaxies could remain unchanged. This could explain 
the observed abundance of large thin bulgeless disks unaffected by major mergers, and in those 
rare mergers between spirals, many tidal dwarf galaxies (TDGs) would be formed. This could 
lead to the possibility that many dwarf galaxies are not primordial but have been formed tidally 
in these encounters (Kroupa et al. 2010). These populations of satellite galaxies, associated with 
globular clusters that formed along with them, would naturally appear in closely related planes, 
thereby providing a natural solution to the satellites phase-space correlation problem (Ibata et 
al. 2013). To explain why these TDG dwarf satellites appear dark matter-dominated, while 
TDGs should not if DM is collisionless (see also problem 2.4.d hereabove), MONDian dynamics 
would be the only possible explanation. What is more, the density-morphology relation for dwarf 
ellipticals could also find a natural explanation (Dabringhausen &: Kroupa 2012). However, all 
this would require a covariant parent theory of MOND, which will also have extreme challenges 
to overcome, and most notably the missing baryon problem that MOND suffers in rich clusters 



of galaxies and the high third peak in the acoustic power spectrum of the CMB. 



4. MOND-inspired covariant theories: carrying their own problems 

We provide hereafter a non-exhaustive list of the most popular covariant theories that have been 
proposed to reproduce the MOND phenomenology on galaxy scales. More details and references 
(as well as a more complete list including proposals based on non-local theories and modified 
inertia) can be found in Famaey & McGaugh (2012), as well as in Bruneton & Esposito-Farese 
(2007) and Skordis & Zlosnik (2012). We also point out hereafter some of the main weaknesses 
of all these approaches. It should be highlighted again that these theories should not be confused 
with the MOND paradigm itself, which is not a theory per se (see e.g. Milgrom 2012). 



4-1. Scalar-tensor k-essence (RAQUAL) 

Main ingredients: An 'Einstein' metric g^ v (with Einstein-Hilbert action), a k-essence scalar 
field (f>, and normal matter from the standard model coupling to a combination of both. 



Action: The matter couples to the 'physical metric': 

9iw = e 2<t> g^ (3) 

and the scalar field is given a so-called k-essence action, with no potential and a non-linear, 
aquadratic, kinetic term: 



c 4 



J d 4 x^f(X), (4) 



where k is a dimensionless constant, / ~ c 2 /ao is a length-scale, and X = kPg^cj),^^^. The 
function f(X) yields the MOND dynamics in the approriate limit X <C 1 (Bekenstein &: Mil- 
grom 1984). 

Main weakness: Does not enhance gravitational lensing in galaxies because of the conformal 
relation between the Einstein and physical metrics, a priori enough to reject the theory. 

4-2. Tensor- Vector-Scalar theory (TeVeS) 

Main ingredients: An Einstein metric g^ v (with Einstein-Hilbert action), a k-essence scalar 
field 4>, a dynamical unit-norm vector field allowing a disformal relation between the Ein- 
stein and physical metrics (thus allowing to enhance gravitational lensing), and normal matter 
coupling to a combination of these fields. 

Action: The matter couples to the physical metric: 

g»» = e- 2,l> g^ - 2amh(2<f>)U„U v , (5) 

the scalar field is given a k-essence action (Eq. and the vector field (Bekenstein 2004, Skordis 
2008): 

Su = -J^q J d A x^g [K^U^U^ - XWUpUv + 1)] , (6) 

where 

K ^nu = Cl gWgP" + c 2 ~g a P~g t * v + c 3 g au g^ + c 4 U a U»f v (7) 



for a set of constants c±, €2,03,04, and A is a Lagrange multiplier enforcing the unit-norm. 



Main weaknesses: (i) Possible stability issues (especially inside matter), (ii) requires a fine- 
tuned form of f(X) to evade strong- field observational constraints in the solar system and binary 
pulsars, (iii) cannot produce a high third peak in the angular power spectrum of the CMB 
without leading to an exaggerated integrated Sachs- Wolfe (ISW) effect, or without resorting to 
additional collisionless hot dark matter (appealing to which might overproduce superclusters 
and massive galaxy clusters as pointed out by Angus & Diaferio 2011), etc. 



4-3. Generalized Einstein- Aether theories (GEA) 

Main ingredients: Normal matter coupling to the physical metric g^ u (with Einstein-Hilbert 
action), and a unit-norm dynamical vector field Uu, also coupling to the physical metric. 

Action: Zlosnik et al. (2006) showed TeVeS to be expressible as a pure Tensor- Vector theory 
in the matter frame, with the physical metric at the same time satisfying the Einstein-Hilbert 
action and coupling minimally to the matter fields, just like in General Relativity (GR). The 
GEA approach makes use of this fact to devise a simpler Tensor- Vector action in the matter 
frame, where the Einstein-Hilbert and matter actions remain precisely as in GR, but with an 
additional unit-norm vector field (Zlosnik et al. 2007): 

c 4 



S u = -T7-?^2 I d x ^9 f(X) - l'\{g» v UpU v + 1) , (8) 
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where (see Eq. [7] and replacing g^ v by g^ v ) 

X = l 2 K a ^Up, a U v>IM (9) 

and 

In the static weak-field limit, the unit-norm constraint fixes the vector field in terms of the 
metric, and from there we have that, in the weak-field limit, X oc | \7<I> | 2 . 



Main weaknesses: Mostly the same as TeVeS for the CMB and ISW (Zuntz et al. 2010) and 
for the additional hot dark matter, BUT, contrary to TeVeS, it can be chosen to approach GR 
as fast as needed for high accelerations (with no fine tuning) so it avoids all the problems TeVeS 
has in binary pulsars, and the unacceptably large preferred- frame effects of TeVeS in the solar 
system. 



44. Bimetric theories (BIMOND) 

Main ingredients: Two metrics g^ and g^ u , normal matter coupling to g^, and a new form 
of 'twin matter' coupling to the second metric g^ u . 

Action: The heart of this class of theories (Milgrom 2009) is an interaction term between the 
two metrics, through a function of scalars constructed from C^ u /ao where = T^ v — f 1 " . For 
instance, one can write: 

c 4 f 

S = 5 m [matter, ff^]+5 m [twin matter, ff^]+ / d 4 x[ay/^g~R + fiJ^gR- 2(gg) 1, H 2 f(X)], 

(11) 

where / = c 2 /ao, and 

v ;2 „u,v fs _ 



X = ?g» v {ClA - C«C$ a ). (12) 



Judiciously chosing the a and j3 parameters yields a whole class of BIMOND theories with 
various phenomenological behaviors. For instance, in matter-twin matter symmetric versions 
(a = j3 = 1), and within a fully symmetric matter-twin matter system, a cosmological constant 
is given by the zero-point of the function /, naturally of the order of 1, thereby naturally lead- 
ing to A ~ do for the large-scale Universe. Matter and twin matter would not interact at all in 
the high-acceleration regime, and would repel each other in the MOND regime (i.e., when the 
acceleration difference of the two sectors is small compared to oq), thereby possibly playing a 
crucial role in the Universe expansion and structure formation. 

Main weaknesses: (i) Should still be checked against the existence of ghost modes, (ii) 
structure formation etc. to be studied in more details, (iii) not clear how these theories can 
explain the angular power spectrum of the CMB (even when appealing to twin matter). 

4-5. Dipolar dark matter (DDM) 

Main ingredients: Normal matter coupling to the GRphysical metric g^ u , dark matter also 
coupling to the physical metric, but carrying a space-likqj four-vector gravitational dipole mo- 
ment 

Action: The dipolar dark matter medium (Blanchet 2007, Blanchet & Le Tiec 2009) is described 
as a fluid with mass current = pu^ endowed with the dipole moment vector £^ (which will 
affect the total density in addition to the above mass density p), with the following action: 



where P is the norm of the projection perpendicular to the four-velocity of the polarization field 
P^ = p£V. The Poisson equation in the weak- field limit is then recovered as: 



From there, in order to reproduce the MOND phenomenology in galaxies, one appeals to a 
'weak-clustering hypothesis', namely the fact that, in galaxies, the dark matter fluid does not 
cluster much (p <C pb) and is essentially at rest (v = 0) because the internal force of the fluid 
precisely balances the gravitational force, in such a way that the polarization field P is precisely 
aligned with the gravitational one g, and g oc —W'{P). At the cosmological level, the monopolar 
density of the dipolar atoms p will play the role of CDM, while the minimum of the potential 
W(P) naturally adds a cosmological constant term A ~ Oq, thus explaining this numerical coin- 
cidence and making the theory precisely equivalent to the ACDM model at linear order for the 
expansion, for large scale structure formation, and for the CMB (naturally explaining the high 
third peak). 

Main weaknesses: (i) Not clear that the weak-clustering hypothesis in galaxies would be 
the only natural outcome of structure formation within this model, (ii) the weak clustering 
hypothesis in itself might be problematic for explaining the residual missing mass in galaxy 
clusters, due to the fact that this residual mass should essentially be concentrated in the central 
parts of these objects, thus leaving the model to rely on baryonic dark matter or additional hot 
dark matter to explain galaxy clusters. 

1 This is to be contrasted with the time-like nature of TeVeS and GEA vector fields in the static weak-field limit 




(13) 



V.(g-47rP) = ^G{ Pb + p). 



(14) 



5. Conclusion 

If one is familiar with cosmology and large scale structure, it must seem rather peculiar that 
anyone would seriously consider alternatives to ACDM based on the MOND phenomenology 
But if one is more concerned with the observed phenomenology in a wide range of galaxy data, 
it seems just as odd to invoke non-baryonic and collisionless CDM together with fine-tuned 
feedback to explain the appearance of an effective force law that appears to act with only the 
observed baryons as a source in galaxies (see Sect 2.4 and 3, see also Famaey & McGaugh 2012 for 
an extensive review). Accepting this, and building a theory to account for this phenomenology, 
could perhaps help resolving (see Sect. 3) many of the other challenges (Sects. 2.1-2.3) currently 
faced by ACDM. Nevertheless, consistent covariant theories that have currently been proposed 
along these lines are at best effective, and bring with them their own challenges (see Sect. 4). 
Indeed, the most important aspect before one rejects any model is to have a 'simpler' model at 
hand, that still reproduces the successes of the earlier favored model but also naturally predicts 
the discrepant data. While MOND, as a paradigm, has a priori predicted a lot of observations 
which ACDM cannot explain in galaxies, it is absolutely fair to say that there is currently no 
alternative which does better overall than ACDM, and in favor of which Ockham's razor would 
be. It would however probably be a great mistake to persistently ignore the fine-tuning problems 
for ACDM and the related uncanny successes of the MOND paradigm on galaxy scales, as they 
could very plausibly point at a hypothetical better new theory of the dark sector. It would for 
instance be extremely exciting if one would manage to find a physical connection between the 
dark energy sector and the galactic MOND phenomenology in the weak-field limit through the 
A ~ Oq coincidence. In this sense, ideas based on entropic gravity are an interesting line of 
thought (Verlinde 2011, Ho et al. 2010, Klinkhamer Sz Kopp 2011, and others, all somewhat 
inspired by Milgrom 1999), but one should nevertheless keep in mind the observational challenges 
that any such new theory of the dark sector could face, especially regarding the third peak of the 
acoustic power spectrum of the CMB. It is, on the other hand, conceivable that a fully successful 
alternative theory does not exist, and that the apparent MONDian behavior of galaxies will be 
explained through small compensatory adjustments of the current ACDM paradigm. But one 
has to realize that this is far from trivial, and one would have to demonstrate how this could 
possibly occur. One also has to realize that most currently proposed MOND-inspired effective 
theories (see Sect. 4) do include new fields, so that the MOND paradigm does not necessarily 
imply the absence of a dark sector. In any case, the existence of a characteristic acceleration 
ao ~ A 1 / 2 playing various roles in many apparently independent galactic scaling relations is 
by now an empirically established fact, and it is thus mandatory for any successful model of 
galaxy formation and evolution to explain it, together with all the other challenges mentioned 
in Sect. 2. 
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